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Abstract 

Integer  linear  programs  arise  in  many  situations,  and  solving  such  problems  can  be 
computationally  demanding.  One  way  to  solve  them  more  efficiently  is  by  exploiting 
the  symmetry  within  their  formulation.  This  paper  proves  that  the  symmetry  group 
for  the  linear  programming  relaxation  of  2-level  orthogonal  array  problems  of  strength 
2  is  a  particular  semi  direct  product. 
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SYMMETRY  GROUPS  FOR  LINEAR  PROGRAMMING 
RELAXATIONS  OF  ORTHOGONAL  ARRAY  PROBLEMS 

I.  INTRODUCTION 


1.1  Motivation 

Solving  integer  linear  programs  (ILPs)  is  a  common  problem,  and  it  is  not  always 
an  easy  task.  In  fact,  doing  so  often  requires  a  substantial  amount  of  compuational 
time.  There  are  several  methods  available  to  solve  ILPs  more  efficiently.  In  particular, 
the  Margot  ILP  solver  is  able  to  decrease  computation  time  by  orders  of  magnitude 
by  taking  advantage  of  symmetry  within  the  formulation  of  a  problem  [4] .  In  order  to 
take  full  advantage  of  the  symmetry  in  a  given  formulation,  one  must  know  as  much 
about  the  symmetry  group  as  possible.  Ideally,  the  entire  symmetry  group  will  be 
known.  Geyer  [2]  developed  an  algorithm  for  finding  the  symmetry  group  for  linear 
programming  (LP)  relaxations  of  ILPs  with  equality  constraints  by  using  projection 
matrices.  The  goal  of  this  research  is  to  prove  a  special  case  of  a  conjecture  resulting 
from  Geyer  [2], 

1.2  Research  Contribution 

This  research  has  proven  that  the  symmetry  group  for  the  LP  relaxation  of  an 
ILP  formulation  of  a  2-level,  k- factor,  strength  2  orthogonal  array  is  A/?  x  S^+i  ■  At 
face  value,  this  result  provides  the  symmetry  group  for  the  LP  relaxation  of  certain 
ILPs.  The  methods  used  to  prove  this  result  may  also  be  useful  for  finding  other 
such  symmetry  groups  with  different  strengths  and  levels.  Furthermore,  this  result 
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provides  a  theoretical  verification  of  Geyer’s  [2]  algorithm.  Finally,  the  methods 
employed  in  this  research  may  be  useful  for  developing  more  efficient  computational 
algorithms. 

1.3  Organization  of  Thesis 

This  thesis  is  divided  into  four  chapters  and  an  appendix.  Chapter  2  contains  a 
brief  discussion  of  some  literature  that  is  directly  related  to  this  research.  Chapter 
3  covers  the  original  work  performed  and  is  divided  into  three  sections.  The  first 
lays  the  groundwork  for  the  work  that  follows.  The  second  section  walks  through  the 
relatively  simple  case  of  a  strength  one  orthogonal  array,  and  the  third  delves  into 
the  strength  two  case.  Chapter  4  outlines  some  topics  for  future  research  that  relate 
to  this  work.  The  appendix  contains  the  code  that  was  used  to  validate  one  of  the 
proofs. 
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II.  LITERATURE  REVIEW 


Experiments  are  generally  intended  to  allow  insight  into  situations.  How  an  ex¬ 
periment  is  designed  is  of  great  importance  to  the  effectiveness  of  the  experiment  and 
the  relevance  of  the  results.  A  factorial  design  is  a  collection  of  factors  that  assume 
a  finite  number  of  level  combinations.  If  one  design  can  be  obtained  from  another  by 
permuting  runs,  factors,  or  levels,  the  two  designs  are  said  to  be  isomorphic.  For  cer¬ 
tain  linear  models,  orthogonal  arrays  are  the  most  efficient  class  of  factorial  designs. 
An  orthogonal  array  with  N  runs  (rows),  k  factors  (columns),  s  levels,  and  strength 
t  is  denoted  OA(iV,  k,  s,  t )  .  The  index  of  such  an  orthogonal  array,  A  ,  is  the  number 
of  times  every  f-tuple  appears  within  each  combination  of  t  columns.  The  index  is 
typically  omitted  from  notation  because  A  =  N/s *  . 

If  an  1LP  contains  any  variables  that  can  be  permuted  without  changing  the 
feasibility  and  optimality  of  its  solutions,  it  is  said  to  be  symmetric.  Margot  [4] 
defined  the  symmetry  group,  G  ,  of  an  1LP  to  be 

G  =  {7T  e  Sn |cTx  =  cr7r(x)  and  7r(x)  eTVxe  IF} 

where  crx  is  the  objective  function  of  the  ILP,  and  J-  is  the  set  of  all  feasible  solutions. 
Symmetric  ILPs  can  arise  from  a  variety  of  problem  formulations.  In  particular,  ILPs 
for  enumerating  orthogonal  arrays  arc  highly  symmetric. 

Optimal  solutions  to  ILPs  are  commonly  found  with  branch- and-bound  or  branch- 
and-cut  algorithms.  In  the  case  of  symmetric  ILPs,  many  of  the  subproblems  in  the 
enumeration  tree  are  isomorphic.  As  a  result,  a  considerable  amount  of  computational 
time  is  wasted  on  solving  identical  problems  repeatedly.  Thankfully,  Margot  [4]  de¬ 
veloped  a  solver  that  is  able  to  decrease  and  potentially  eliminate  such  redundant 
computations  by  exploiting  a  subgroup  of  an  ILP’s  symmetry  group  when  pruning 
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its  enumeration  tree.  Exploiting  larger  subgroups  results  in  increased  reductions  in 
redundant  computations,  and  the  greatest  reduction  is  attained  if  the  full  symmetry 
group  of  the  ILP  is  exploited.  Hence,  it  is  desirable  to  find  larger  subgroups  of  the 
symmetry  group  of  an  ILP. 

For  an  ILP  with  objective  function  c1  x  ,  m  x  n  constraint  matrix  A  ,  and  righthand 
side  b  with  0  <  x  <  d  ,  let  A(7t,  a)  be  the  matrix  found  by  permuting  the  rows  and 
columns  of  A  according  a  and  tt  ,  respectively.  The  automorphism  group  of  such  an 

ILP  is 

G(A,  b,  c,  d)  =  {n  e  ^^(c)  =  c,  7r(d)  =  d,  and  €  Sm  :  A(i r,  a)  =  A,  cr(b)  =  b} 

where  a  preserves  equalities  and  inequalities.  Clearly,  G( A,  b,  c,  d)  <  G  ,  and  com¬ 
putational  experiments  suggest  that  efficiency  can  be  improved  by  several  orders  of 
magnitude  if  Margot’s  solver  is  used  with  G( A,  b,  c,  d)  on  highly  symmetric  ILPs  [4], 
Because  the  makeup  of  the  symmetry  group,  G  ,  is  determined  by  the  ILP’s 
feasible  set,  identifying  all  symmetries  in  an  ILP  is  quite  difficult  and  remains  an 
open  problem.  Margot  [4]  proved  that  deciding  if  G  =  Sn  is  NP-Complete.  Therefore, 
finding  G  for  any  given  ILP  is  NP-Hard.  In  order  to  find  many  of  the  symmetries  in 
an  ILP,  one  can  simply  find  the  symmetry  group  of  the  LP  relaxation  of  the  ILP.  The 
LP  relaxation  symmetry  group,  GLP  is  the  set  of  all  permutations  of  variables  that 
send  LP  feasible  points  to  LP  feasible  points  with  the  same  objective  function  value. 
For  an  LP  relaxation  without  equality  constraints  where  each  constraint  in  Ax  <  b 
is  a  facet  (non-redundant),  GLP  =  G(A,b,c,d)  .  Let  iF(k,s,t )  and  Jr(/c,s,f)LP  be 
the  sets  of  feasible  solutions  of  an  ILP  and  its  LP  relaxation,  respectfully.  Note  that 
C  F(k,s,t)LP  ,  so  GLP  <  G  . 

For  an  LP  relaxation  with  equality  constraints,  G( A,  b,  c,  d)  C  GLP  .  Geyer  [2] 
developed  an  algorithm  for  finding  GLP  of  such  an  LP.  Furthermore,  Geyer  [2]  ob- 
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served  that  G(A,b,c,d)  C  GLP  for  the  LP  relaxation  of  an  ILP  formulation  with 
equality  constraints  for  finding  OA(N,k,2,t)  when  t  is  even.  Geyer,  Bulutoglu,  and 
Rosenberg  [3]  explicitly  found  G(A,  b,  c,  d)  for  this  formulation.  They  also  found  GLP 
for  a  formulation  without  equality  constraints.  This  thesis  verifies  Geyer’s  compu¬ 
tational  observations  by  explicitly  finding  GLP  for  an  ILP  formulation  with  equality 
constraints.  The  equality  constraints  of  this  ILP  formulation  are  linear  combinations 
of  those  of  the  ILP  formulation  with  equalities  in  [3].  Furthermore,  both  ILPs  have 
the  same  number  of  non-redundant  equality  constraints.  Hence,  one  can  go  back  and 
forth  between  the  two  ILP  formulations  by  applying  a  sequence  of  row  operations  to 
the  equality  constraints  of  each.  This  implies  that  the  feasible  sets  of  the  LP  relax¬ 
ations  of  these  ILP  formulations  are  the  same,  so  GLP  must  also  be  the  same  for  the 
two  ILPs. 

The  ILP  formulation  used  in  this  thesis  stems  from  the  concept  of  J-characteristics. 
Let  the  frequency  vector,  f  ,  of  a  2-level  factorial  design,  D  ,  have  the  frequency  of 
each  of  the  2k  possible  factor  level  combinations  as  its  entries.  Hence,  f  determines  D 
up  to  reordering  of  factor  level  combinations.  For  a  2-level  design,  D  ,  with  N  runs 
(factor  level  combinations)  and  k  factors,  J-characteristics  are  given  by 

N 

Ji  =  dij 

*= i  jez 

for  l  C  Z*,  .  It  has  been  shown  that  D  is  uniquely  determined  by  its  J-characteristics 
up  to  reordering  of  its  runs;  furthermore,  D  is  an  orthogonal  array  of  strength  t  if 
and  only  if  Ji  =  0  for  all  l  C  7Lk  with  |/|  <t  ,  where  l  ^  0  [7], 
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III.  ORIGINAL  WORK 


3.1  Preliminary  Steps 

Let  1  be  the  column  vector  of  length  2k  for  which  every  entry  is  one.  For  i  = 
1, . . .  ,k  ,  let  Xj  be  the  ith  column  of  the  /e-factor,  2-level  (±1)  full  factorial  design. 
For  distinct  G  {1, . . . ,  k}  ,  let  Xj^...^.  represent  the  j-factor  interaction  term 

given  by  the  Hadamard  product  Xjx  o  •  •  •  o  x* .  . 

Consider  the  equation 

Mf  =  J  (1) 

t 

from  [7],  where  M  is  the  ^  (*?)  by  2k  matrix 

i= 0  1 


1T 


X 


T 

1 


M  = 


x 


T 

k 


J  is  the  quasi-lexicographically  ordered  J-characteristic  vector  with  entries  J/  for 
|/|  <  t  ,  and  f  is  the  frequency  vector  of  a  hypothetical  OA(N,k,2,t)  .  Then  by  the 
result  in  [7] 


N 

0 


0 
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Our  goal  is  to  find  the  subgroup  of  the  permutation  group  S2k  that  sends  feasible 
solutions  (f  G  Q>0)  to  feasible  solutions.  This  is  the  symmetry  group  of  the  LP 
relaxation  of  the  ILP  based  on  Equation  1,  where  the  objective  function  is  taken  to 
be  the  zero  vector.  The  equality  constraints  of  this  ILP  are  linear  combinations  of 
those  of  the  orthogonal  array  defining  ILP  in  [3] .  Both  ILPs  have  the  same  inequality 

t. 

constraints,  and  each  ILP  has  0)  non- redundant  equality  constraints.  Hence,  both 

i= 0 

ILPs  have  the  same  LP  relaxation  feasible  set,  and  this  implies  that  both  have  the 
same  LP  relaxation  symmetry  group.  From  this  point  on,  we  shall  refer  to  this  group 
as  G  . 

Theorem  1  The  symmetry  group  G  is  precisely  the  intersection  of  the  automorphism 
group  of  the  row  space  o/M  and  the  permutation  group  S2k  ,  written  Aut(Row(M))  D 
S2k  .  That  is,  G  is  the  set  of  permutations  that  preserve  Row( M)  . 

Proof  Observe  that 

N_ 

2  k 

p  _ 

N_ 

2k 

is  a  particular  solution  to  Equation  1.  As  such,  every  solution  f  can  be  written  in  the 
form  f*  +  f'  where  f'  G  NullfWl)  .  Let  g  G  G  be  arbitrary.  Then  g(f)  is  a  solution  to 
Equation  1.  That  is, 

Mg(i)  =  M  g({*  +  f')  =  J 

Because  g  G  G  <  S2k  , 

M[<?(f*)  +  g(f')\  =  M[f*  +  g(f7)]  =  J 

and  thus 

MP  +  =  J  +  Mg(f')  =  J 
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Therefore, 


M(/(f/)  =  0 

so  we  see  that  g( f')  G  Null(Wl)  which  means  g  must  preserve  Null(WL)  .  Because 
g  G  G  was  arbitrary,  we  have  shown  G  <  Aut(Null( M))  D  S2k  . 

Now  let  h  G  Aut(Null( M))  D  S2k  be  arbitrary.  Then 

M/i(f)  =  M/i(f*  +  f') 

=  M[/i(f*)  +  h{i')} 

=  M[f*  +  /i(f')] 

=  Mf  +  M/i(f') 

=  J  +  0 
=  J 

Hence,  h  G  G  ,  and  because  h  was  arbitrary,  Aut(Null( M))  D  S2k  <  G  .  Noting 
that  Aut(Null(Wl))  =  Aut(Row(WL ))  ,  we  conclude  that  G  =  Aut(Null(Wl ))  D  S2k  = 
Aut(Row  (M) )  fi  S2k  . 

3.2  The  Symmetry  Group  for  the  Strength  One  Case 

Before  investigating  the  symmetry  group  for  a  strength  two  orthogonal  array,  it 
is  only  natural  that  we  should  address  the  strength  one  case.  As  such,  for  now  we 
assume  an  OA(iV,  k,  2, 1)  in  Equation  1.  Let 

B  =  { 1 ,  x1; . .  •  ,xfc} 

Clearly,  B  is  an  orthogonal  basis  for  ifotu(M)  .  For  all  g  E  G  ,  we  know  g(B)  must  also 
be  an  orthogonal  basis  for  ifotu(M)  because  by  Theorem  1,  g  G  Aut(Row(WL))  D  S2k  , 
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and  elements  of  S2k  preserve  angles.  Furthermore,  for  every  x  e  B  ,  we  know  that 
g(x)  can  be  represented  uniquely  as  a  linear  combination  of  the  elements  of  B  .  That 
is, 

g(x)  =  Aol  +  AiXi  +  •  •  •  +  AfcXfc  (2) 

Lemma  1  Let  x  e  B  .  If  x  —  1  in  Equation  2,  then  Ao  =  1  ,  and  \  =  0  for 
i  —  1, . . . ,  k  .  Otherwise,  A0  =  0  . 

Proof  Suppose  x  =  1  .  Because  g  G  G  <  S2k  ,  g(  1)  =  1  which  uniquely  satisfies 
Equation  2.  For  i  —  1, . . . ,  k  ,  we  know  g(xf)  must  be  orthogonal  to  g(  1)  =  1  ,  so 
Ao  =  0  whenever  x^  1  . 

Lemma  2  If  {x', , . . . ,  x'k }  is  obtained  from  {xi, . . .  ,  x*,}  ,  the  columns  of  the  full 
factorial  2k  design,  by  permuting  rows,  then  there  exists  a  permutation  cr  e  such 
that  for  all  i  e  {1, . . , ,  k}  satisfying  x'  e  span{x i, . . . ,  x^,)  ;  x'  =  ±xCT(q  . 

Proof  Let  {x'j , . . . ,  x'k\  is  obtained  from  {x1? . . .  ,x*.}  ,  the  columns  of  the  full  fac¬ 
torial  2k  design,  by  permuting  rows.  Suppose  x'  e  span(xi, . . . ,  xk)  .  Then  x'  = 
AiXi  +  •  •  •  +  AfcXfc  ,  and  we  have  the  system  of  equations 

Ai  +  ■  ■  ■  +A  k  —  il 

Ai  +  •  •  •  —A  k  =  ±1 

— A|  —  •  •  •  — Afc  =  ±1 

Subtracting  the  second  equation  from  the  first  equation  gives  A*,  e  {0,  ±1}  .  Choosing 
other  pairs  of  equations  similarly  yields  A j  G  {0,  ±1}  for  j  —  1, , . . ,  k  .  Because  B  is 
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an  orthogonal  set,  the  Pythagorean  Theorem  gives 

Nil2  =  EllV.ll2 
.=1 
k 

=  £(Miwi)2 

.=i 

=  DAPIMI2 
.'= 1 

=  X>?IMI2 
.'= 1 

Because  row  permutations  are  norm-preserving,  1 1 x ■  |  =  2k  =  ||xj||  for  j  =  1 , ,k  . 
Thus, 

k 

3= 1 

Because  Aj  <  {0,  ±1}  for  j  =  l,...,/c  ,  there  is  exactly  one  nonzero  A  j  G  {±1}  , 
and  x-  =  ±Xj  .  Row  permutations  also  preserve  orthogonality,  so  for  every  distinct 

i  G  { 1 _ ,  }  such  that  x'  G  span(x i, . . .  ,x^)  ,  there  is  a  unique  j  G  {1, . . . ,  k} 

satisfying  x-  =  ±Xj  .  Thus,  there  exists  a  permutation  a  G  Sk  such  that  x-  =  ±xcr(q 
for  all  i  G  {1, . . . ,  k}  such  that  x'  G  span(x.i , . . . ,  xfe)  . 

Lemma  3 

|G|  <  2 kk\ 

Proof  Let  g  G  G  be  arbitrary.  From  Lemma  1,  we  have  that  g(  1)  =  1  ,  and 
gf(xj)  =  AiXi  +  •  •  •  +  AfcXfc  for  i  =  1, . . . ,  k  .  Now  by  Lemma  2  we  know  there  exists  a 
permutation  cr  G  Sk  such  that  x'  =  ±xCT(j)  for  i  —  1, . . . ,  k  .  That  is,  g  is  essentially  a 
signed  permutation  of  the  k  main  effects,  so  g  is  one  of  at  most  2 kk\  elements  in  G  . 

Next,  by  finding  a  subgroup  of  G  that  attains  the  upper  bound  on  size,  we  de¬ 
termine  the  size  and  the  structure  of  G  .  We  first  introduce  some  terminology  and 
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notation  as  given  by  Rotman  [6]  on  pages  167  and  172. 


Definition  Let  K  be  a  (not  necessarily  normal)  subgroup  of  a  group  G  .  Then  a 
subgroup  Q  <  G  is  a  complement  of  K  in  G  if  K  fl  Q  —  1  and  KQ  =  G  . 

Definition  A  group  G  is  a  semidirect  product  of  K  by  Q  ,  denoted  by  G  —  K  x  Q  , 
if  K  <  G  and  K  has  a  complement  O'  =  0  .  One  also  says  that  G  splits  over  K  . 

Definition  Let  D  and  Q  be  groups,  let  O  be  a  finite  Q-set,  and  let  K  =  ILsi  A,  , 
where  Du  =  D  for  all  u  6  H  .  Then  the  wreath  product  of  D  by  Q  ,  denoted  by 
D  \  Q  (or  by  D  wr  Q),  is  the  semidirect  product  of  K  by  Q  ,  where  Q  acts  on  K  by 
q(du)  =  (dqu)  for  q  G  Q  and  (du)  G  n^eo  Aj  •  The  normal  subgroup  K  of  D  l  Q  is 
called  the  base  of  the  wreath  product  . 

In  the  factorial  desing  setting,  Sj.  is  the  permutation  group  of  k  factors.  The 
multiplicative  group  {±1}  that  multiplies  columns  is  isomorphic  to  S2  ■  Naturally, 
Sf  is  the  direct  product  of  k  copies  of  S2  ■  We  now  see  that  Sf  is  the  base  of  S2  l  Sk  . 

Lemma  4 

sk2  <  S2 1  sk 

Proof  Per  the  definition  of  the  wreath  product,  acts  on  Sk  by  permuting  the 
entries  of  each  {±l}fc  vector  in  S2  ■  Let  (p  G  S2  and  a  G  Sk  be  arbitrary.  Clearly, 
a  (pa G  S2  .  Furthermore,  a~x(pa  G  S2  ,  and  a(a~1(pa)a~1  =  (p  .  Thus,  aS2a~l  =  S2 
for  all  a  E  Sk  ,  so  S2  <  S2  x  Sk  =  S2l  Sk  . 

This  wreath  product  is  the  set  of  all  signed  permutations  of  x*  for  i  =  1 , ,k 
from  the  full  factorial  2k  design,  where  xj1, . . . ,  xjP  constitute  rows  of  M  .  We  shall 
see  that  this  group  is  a  subgroup  of  G  .  Hence,  by  Lemma  3,  it  is  G  . 
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Theorem  2 


G  =  S2l  Sk 


Proof  Consider  an  arbitrary  element  of  S2l  Sk  .  We  know  it  can  be  written  in  the 
form  0a  where  0  G  Sk  and  a  G  Sk  thanks  to  Lemma  4.  Clearly,  permuting  the  k 
rows  xf , . . . ,  x^  of  M  will  preserve  the  full  factorial  2k  design  as  will  negating  any 
of  these  k  rows.  Hence,  0a  G  S2k  .  Furthermore,  the  signed  permutation  0a  clearly 
preserves  -Row(M)  ,  so  0a  G  Aut(Row( M))  .  Because  0a  G  S2  l  Sk  was  arbitrary,  we 
know  S2l  Sk  <  G  =  Aut(Row( M))  D  S2k  .  Finally,  \S2  l  <Sfc|  =  |S'f||S'/c|  =  2 kk\  which 
is  the  upper  bound  for  \G\  ,  so  G  must  be  exactly  S2l  Sk  ■ 


Corollary  1 


|G|  =  2  kk\ 


Proof  As  a  direct  result  of  Theorem  2,  we  have  \G\  =  S2  l  50  =  2 kk\  . 


3.3  The  Symmetry  Group  for  the  Strength  Two  Case 

From  this  point  forward,  we  assume  an  OA(iV,  k ,  2,  2)  in  Equation  1.  As  in  the 
strength  one  case,  let 


B  =  {1,  X1?  .  .  .  ,  Xfc,  Xi;2,  •  •  • ,  xfc_liA.} 

Once  again,  B  is  an  orthogonal  basis  for  .Roic(M)  ,  and  for  all  g  £  G  ,  we  know  g(B ) 
must  also  be  an  orthogonal  basis  for  i?ou;(M)  because  g  G  Aut(Row( M))  D  S2k  as 
given  by  Theorem  1.  Also,  for  every  x  G  B  ,  we  know  that  g(x)  can  be  represented 
uniquely  as  a  linear  combination  of  the  elements  of  B  .  In  this  case, 


g(x)  —  AqI  +  AiX!  +  •  •  •  +  AfcXfc  +  Ai^xx^  +  •  •  •  +  Xk-i,k^-k-i,k  (3) 
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Similar  to  the  strength  one  case  above,  we  will  arrive  at  the  conclusion  that  for 
any  x65,  every  A  in  Equation  3  must  be  zero  except  for  one,  which  must  have  an 
absolute  value  of  one.  The  following  several  lemmas  serve  to  lead  us  to  this  conclusion. 

Lemma  5  Let  x  G  B  .  If  x  —  1  in  Equation  3,  then  Ao  =  1  ,  and  A*  =  0  for 
i  —  1, . . . ,  k,  (1, 2), . . . ,  (k  —  1,  k)  .  Otherwise,  A0  =  0  . 

Proof  Suppose  x  =  1  .  Because  g  G  G  <  S2k  ,  g(  1)  =  1  which  uniquely  satisfies 
Equation  3.  For  i  —  1, . . . ,  k,  (1,  2), . . . ,  (k  —  1,  k)  ,  we  know  g(xj)  must  be  orthogonal 
to  g(  1)  =  1  ,  so  Ao  =  0  whenever  x/1  . 

Lemma  6  Let  x  G  B  .  If  x  ^  1  in  Equation  3,  then  \  G  {0,  ±0.5,  ±1}  for  i  = 
1,...,*, (1,2),-.. ,(*-!,*)  ■ 

Proof  Suppose  x  ^  1  .  Then  Equation  3  becomes  g{x)  =  AiXi  +  •  •  •  ±  X^x-k  ± 
Ai,2xi,2  +  •  •  •  ±  Xk-i}k^k-i,k  ■  Because  these  basis  vectors  are  the  columns  of  the  full 
factorial  2k  design  and  the  corresponding  2-factor  interactions  obtained  by  taking  the 
appropriate  pairwise  Hadamard  products  of  the  individual  columns  (main  effects),  we 
have  the  system  of  equations 


Ai 

±  •• 

■  •  ±Afc 

±Al,2 

±  •• 

■  •  ±Aiy- 

±  •• 

■  •  +Xk~i,k  — 

±1 

Ai 

_  . . 

A  k 

—  Al,2 

_  .  . 

■  ■  —  Ai,*; 

±  •• 

■  •  ±Afc_i,fc  = 

±1 

-Ai 

±  •• 

'  •  ±Afc 

—  Al,2 

—  •  • 

'  •  — Ai^ 

±  •• 

■  •  ±Afc-i,jt  = 

±1 

-Ai 

■  •  —Xk 

±Al,2 

±  •• 

■  •  ±Ai,fc 

±  •• 

■  •  ±Afc-i,fc  = 

±1 

Subtracting  the  last  equation  from  the  first  gives  AH - b  A*,  G  {0,  ±1}  .  Taking  the 

difference  of  the  middle  equations  likewise  provides  Ai - A k  G  {0,  ±1}  .  Summing 

these  two  expressions  results  in  the  conclusion  Ai  G  {0,  ±0.5,  ±1}  .  Choosing  other 
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sets  of  equations  similarly  yields  \  G  {0,  ±0.5,  ±1}  for  i  =  1, . . . ,  k,  (1,  2), . . . ,  (k  — 

M')  • 

Lemma  7  Let  x  G  B  and  g  G  G  .  If  x  ^  1  in  Equation  3,  then  either  g(x)  =  ±x* 
for  i  G  {1, . . . ,  k,  (1, 2)|, . . ,  (k  —  1,  k)}  or  g(x)  =  ±0.5xa  ±  0.5x(,  ±  0.5xc  ±  0.5xrf  for 
distinct  a,b,c,d  G  {1, . . . ,  k,  (1,  2), . . . ,  (k  —  1,  A;)}  . 

Proof  Suppose  x^  1  .  Because  B  is  an  orthogonal  set,  the  Pythagorean  Theorem 
gives 

n<?(x)ir  =  ^iia^h2 

i 

i 

=  iA*i2iixiii2 

i 

=  £A?IMI2 

i 

We  note  that  g  G  S2k  is  norm-preserving,  so  ||g(x)||  =  2k  =  ||x,;||  for  i  = 

1, . . . ,  &,  (1,  2), . . . ,  (k  —  1,  k)  .  Thus, 

EA±! 

i 

Clearly,  not  every  A*  can  be  zero.  If  A i  G  {±1}  for  some  i  G  {1, . . . ,  A;,  (1,  2 (k  — 
1,  k)}  ,  then  \j  =  0  for  all  j  G  {1, . . . ,  k,  (1,  2), . . . ,  (k  —  1,  k)}  such  that  i  ^  j  . 
Otherwise,  there  must  be  distinct  a,b,c,dE  {1, . . . ,  k,  (1, 2), . . . ,  (k  —  1,  k )}  such  that 
Aa,  A;>,  Ac,  A d  G  {±0.5}  ,  and  every  other  A  is  zero.  That  is,  either  g(x)  =  ±x,:  for 
%  G  {1, . . . ,  k,  (1,  2), . . . ,  (k  —  1,  k )}  or  g(x)  =  ±0.5xa±0.5x5±0.5xc±0.5xrf  for  distinct 
a,  h,  c,  d  G  {1, . . . ,  k,  (1,  2), . . . ,  (k  -  1,  A;)}  . 

Lemma  8  Let  x  G  B  and  g  G  G  .  If  r/(x)  is  of  the  second  form  given  in  Lemma  7, 
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then  g(x)  =  ±0.5xaife  ±  0.5xaiC  ±  0.5xb  ±  0.5xc  for  distinct  a,  b,  c  G  {1, . . . ,  k}  . 

Proof  Suppose  g(x)  =  ±0.5xa  ±  0.5x5  ±  0.5xc  ±  0.5xd  for  distinct  a,  b,  c,  d  G 
{1, . . . ,  k,  (1,  2), . . . ,  {k— 1,  k)}  .  Clearly,  xa  ,  xb ,  xc ,  and  xd  cannot  all  be  main  effects, 
for  the  full  factorial  design  will  ensure  some  entry  of  g(x)  equals  2  ^  {±1}  •  Therefore, 
at  least  one  2-factor  interaction  must  be  present  in  the  linear  combination.  Because 
there  are  more  such  linear  combinations  than  would  be  prudent  to  check  manually,  we 
take  advantage  of  R  software  [5]  at  this  stage  of  the  proof.  The  code  used  for  this  step 
is  contained  in  the  appendix.  By  creating  every  essentially  unique  linear  combination 
containing  at  least  one  2-factor  interaction  term  and  checking  whether  they  satisfy 
a  basic  requirement,  we  rule  out  all  possibilities  except  those  of  one  particular  form. 
Specifically,  by  ruling  out  each  linear  combination  where  the  minimum  and  maximum 
entries  in  the  resulting  vector  are  not  —1  and  1  ,  respectively,  we  eliminate  all  linear 
combinations  except  those  of  the  form  ±0.5xfli5  ±  0.5xajC  ±  0.5x5  ±  0.5xc  for  distinct 
a,  b,  c  G  {1, . . . ,  k}  . 

It  is  clear  that  k  >  3  in  order  for  the  form  in  Lemma  8  to  be  viable. 

Lemma  9  Let  x  G  B  and  g  G  G  .  If  for  every  i  G  {1, . . . ,  k}  ,  g(x,:)  ^  ±0.5xa5  ± 
0.5xa  c  ±  0.5x5  ±  0.5xc  for  distinct  a,  5,  c  G  {1, . . .  7  k}  ,  then  g(x)  cannot  be  of  the 
form  in  Lemma  8. 

Proof  Recall  from  Lemma  5  that  g(  1)  =  1  .  Suppose  that  for  every  i  G  {1, . . . ,  k}  , 
we  have  g(x.i)  ^  ±0.5xa  5  ±  0.5xaiC  ±  0.5x5  ±  0.5xc  for  distinct  a,  5,  c  G  {1, . . . ,  k}  . 
Then  we  know  from  Lemmas  7  and  8  that  for  every  i  G  {1,. ,,,  k}  ,  there  exists 
some  j  G  {1, . . . ,  k,  (1,  2), . . . ,  (k  —  1,  &)}  such  that  g(x* )  =  ±Xj  .  Because  g  preserves 
Hadamard  products,  for  every  i  G  {(1,2),...,  (A;  —  1  ,k)}  ,  there  exists  some  j  G 
{1, . . . ,  k,  (1,  2), . . . ,  {k  —  1,  k)}  such  that  g(x,;)  =  ±Xj  .  Hence,  for  every  x  G  B  , 
g(x)  ^  ±0.5xfli5  ±  0.5xajC  ±  0.5x5  ±  0.5xc  for  distinct  a,  b,  c  G  ,k}  . 
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Lemma  10  Let  g  E  G  .  If  for  some  i  E  {1, . . . ,  k}  ,  gigti)  =  ±0.5xa  fe  ±  0.5xajC  ± 
0.5x;,±0.5xc  for  distinct  a,  b,  c  G  {1, . . . ,  k}  ,  then  there  must  exist  some  j  E  {1, . . . ,  k} 
with  i  ^  j  such  that  g(x.j)  =  ±0. 5xa/^  ±0.5xa/)C/  ±0.5x5/  ±0.5xc/  for  distinct  a',  b' ,  c'  G 
{1  ,■••,&}  • 

Proof  Suppose  there  exists  some  i  G  {1, . . . ,  k}  ,  g(x,;)  =  ±0.5xaj5±0.5xajC±  0.5x5  ± 
0.5xc  for  distinct  a,  6,  c  G  {1, . . . ,  k}  .  By  way  of  contradiction,  suppose  there  is  no 
j  G  {1 ,k}  with  i  ^  j  such  that  g(xj)  =  iO-Sx^/  ±  0.5xa/jC/  ±  0.5x5/  ±  0.5xc/ 
for  distinct  a',b',c'  G  {l,...,/c}  .  Then  from  Lemma  6  we  know  that  for  every 
j  E  {1, ...  ,k}  with  i  j  there  exists  some  l  E  {1, . . . ,  k,  (1,  2), . . . ,  (k  —  1,  fc)}  such 
that  gigtj)  =  ±Xi  .  Because  g  preserves  Hadamard  products,  g(x,;)  o  ±x;  must  also 
take  on  a  viable  form,  and  this  implies  that  x^  G  {xa,X5iC}  .  There  can  only  be  one 
such  Xj  because  if  there  were  more  than  one,  their  Hadamard  product  would  be  sent 
to  something  in  {±l,±xa)5iC}  .  But  that  means  only  two  main  effects  (xj  and  Xj) 
get  sent  to  viable  forms,  which  contradicts  Lemma  7.  Thus,  there  must  exist  some 
j  E  {1, ...  ,k}  with  i  jtz  j  such  that  gigtj)  =  TO-Sx^/  ±  0.5xa/jC/  ±  0.5x5/  ±  0.5xc/  for 
distinct  a',  6',  d  E  {1, . . . ,  k}  . 

Lemma  11  Let  g  E  G  .  If  there  exist  distinct  i,j  E  {1, . . . ,  k}  such  that  g(x,:)  = 
±0.5xai5±0.5xajC±0.5x5±0.5xc  for  distinct  a,b,c  E  {1, . . . ,  k}  and  g(x.j)  =  TO-Sx^/i 
0.5xa/iC/  ±  0.5x5/  ±  0.5xc/  for  distinct  a',  br,  d  E  {1, . . . ,  k}  #  then  {a,  b,  c}  =  {a',  b',  c'}  . 

Proof  Suppose  there  exist  distinct  i,j  E  {1  ,...,k}  such  that  g(x*)  =  ±0.5xa)5  ± 
0.5xaiC±0.5x5±0.5xc  for  distinct  a,b,c  E  {1, . . . ,  k}  and  g(x.j)  =  TO.Sxa/^/iO.Sxa/^/i 
0.5x5/  ±  0.5xc/  for  distinct  a b',  d  E  {1, ...  ,k}  .  We  proceed  by  way  of  contradiction 
and  suppose  that  {a,  6,  c}  ^  {afbfd}  .  That  is,  \{a,b,c}  D  {a',b',d}\  <  3  .  We 
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observe  that 


g{xij)  —  ±0.25x0^^/^  ±0.25xai6)(l/jC,  ±0. 25x^6/  ±0.25xa){,c/ 
±0.25xac>a/ib/  ±0.25xa)C)a/)C/  ±0.25xOiC)6/  ±0.25xacc/ 
±0.25x6^/^/  ±0.25xf,ia/!C/  ±0.25x6^/  ±0.25x^c/ 

±0.25  ~3tc,a',v  ±0.25xCia/jC/  ±0.25xCife/  ±0.25xCiC/ 


(4) 


(Case  1:  |{a,  b,  c}  fl  {a',  b',  c'}|  =  0)  Equation  4  clearly  is  not  of  a  valid  form. 
(Case  2:  \{a,b,  c}  fl  {a/,6/,c/}|  =  1)  If  a  ^  a'  ,  4-factor  interaction  terms  will 
remain  in  Equation  4,  so  it  will  not  be  of  a  valid  form.  Suppose  a  =  a'  .  Even  if  the 
3-factor  interaction  terms  were  to  cancel,  the  remaining  2-factor  interaction  terms  are 
insufficient  for  Equation  4  to  be  of  a  valid  form. 

(Case  3:  |{a,  b,  c}  fl  {a',b',c'}\  —  2)  If  a  ^  a!  ,  at  least  one  4-factor  interaction 
term  will  remain  in  Equation  4,  so  it  will  not  be  of  a  valid  form.  Suppose  a  =  a ’  . 
Without  loss  of  generality,  also  suppose  b  =  b'  .  Even  if  the  3-factor  interaction  terms 
and  the  1  terms  were  to  cancel,  the  remaining  2-factor  interaction  terms  and  main 
effect  terms  are  insufficient  for  Equation  4  to  be  of  a  valid  form. 

This  contradicts  Equation  4  being  of  a  valid  form,  so  we  conclude  that  {a,  b,  c }  = 

{«',  V,  c'}  ■ 

Lemma  12  If  k  >  4  ,  and  if  x  ^  1  in  Equation  3,  then  g(x)  =  ±x,:  for  i  e 
{ 1 , . . . ,  A;,  (1,  2),  ...,(&  —  1,  &)}  . 

Proof  Let  k  >  4  ,  and  x  ^  1  in  Equation  3.  By  way  of  contradiction,  suppose 
g(x)  =  ±0.5xa  b±0.5xa  c±0.5xb±0.5xc  for  distinct  a,  b,  c  G  {1, . . . ,  k}  .  From  Lemma 
9,  we  know  there  must  be  some  i  G  {1, . . . ,  k}  such  that  g(xj)  =  TOffix^/  ±0.5xa/j(y  ± 
0.5xfc/  ±  0.5xc/  for  distinct  a',  b',  d  e  {1, . . , ,  k}  .  Lemma  10  guarantees  there  will  be 
another  main  effect  sent  to  a  similar  form  by  g  ,  and  Lemma  11  tells  us  it  will  be  built 
from  the  same  three  distinct  main  effects  and  their  three  distinct  2-factor  interactions. 
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At  most  one  main  effect  could  be  sent  by  g  to  a  form  other  than  that  just  described 
as  noted  in  the  proof  of  Lemma  10.  Now  we  have  at  least  one  more  main  effect  to 
consider,  and  it  must  be  sent  to  a  form  similar  to  that  given  above  and  also  built  from 
the  same  three  main  effects  and  their  three  distinct  2-factor  interactions.  But  now  we 
have  four  main  effects  that  are  sent  to  linear  combinations  of  six  orthogonal  vectors, 
and  the  six  resulting  2-factor  interactions  will  necessarily  also  be  sent  by  g  to  linear 
combinations  of  those  same  six  orthogonal  vectors  (owing  to  the  properties  of  the 
Hadamard  product).  This  means  that  the  ten  new  vectors  cannot  all  be  orthogonal, 
which  contradicts  g  G  G  .  Hence,  g(x)  ^  ±0.5xafe  ±  0.5xa)C  ±  0.5xb  ±  0.5xc  for 
distinct  a,  b,  c  G  {1,  .  Now  from  Lemmas  8  and  7,  we  have  g(x)  =  ±x*  for 

i  G  {l,...,/c,(l,2),...,(/c-  1  ,k)}  . 

Lemma  13  Let  k  >  4  .  Then 

\G\  <  2k{k  +  1)! 

Proof  Let  g  G  G  be  arbitrary.  Note  that  because  g  preserves  Hadamard  products, 
knowing  how  it  acts  on  the  main  effects  will  determine  how  it  acts  on  all  of  B  .  From 
Lemma  12,  we  know  g{x.i)  =  ±Xj  for  i  G  {1, . . . ,  k,  (1,  2), . . . ,  (k  —  1,  k)}  .  Because 
g(xi)2)  must  be  of  a  similar  form,  the  possibilities  for  g(x 2)  are  restricted  depending 
upon  g(xi)  .  If  g(xi)  =  ±x*  for  i  G  {1, , . . ,  k}  ,  then  g(x2)  =  ±x;  for/G{l,...,i  — 
1,  i  +  1, . . ,  ,  k,  (1,  i),  > . . ,  (i  —  1,  i),  (i,  i  +  1), . . . ,  (i,  k)}  .  Otherwise,  ^(xx)  =  ±Xjj  for 
i  <  j  and  i,j  G  {1, . . . ,  k}  ,  so  g(x2)  =  ±x;  for  l  G  {i,  j,  (1,  j), . . . ,  (i  -  l,i),  (i - 

1  ,j),  (*,*  +  1),  (*  +  1,  j),  ■  ■  ■  y(hj  ~  1),  (J  -  1  ,j),  (i,j  +  1),  (jJ  +  !),•••,  (*,  k),  ( j ,  k)}  . 

To  determine  how  many  distinct  possibilities  exist,  we  shall  consider  four  cases,  based 
on  the  forms  of  ^(xi)  and  g(x2)  ,  respectively. 

(Case  1:  main  effect,  main  effect)  Suppose  g(xi)  =  ±x,:  for  i  G  {l,...,/c}  and 
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g(x2)  =  ±x;  for  l  E  {1, . . . ,  i  —  1,  i  +  1, . . . ,  k}  .  Then  there  are  2k  possibilities  for  g(*-i) 
and  2(k  —  1)  for  g(x2)  .  All  of  the  (k  —  2)  remaining  main  effects  must  be  sent  to  pins  or 
minus  the  other  (k  —  2)  main  effects.  That  is,  there  are  (2k)  (2  (k  —  l))(2k~2(k  —  2)!)  = 
2kk\  distinct  possibilities. 

(Case  2:  main  effect,  2-factor  interaction)  Suppose  g(x i)  =  ±x*  for  i  E  {1, . . . ,  k} 
and  g(x2)  =  ±x^  for  l  €  {(1,  i), . . . ,  (i  —  1,  i),  (i,  i  +  1), . . . ,  (i,  k)}  .  Then  there  are 
2k  possibilities  for  g(x i)  and  2 (k  —  1)  for  g(x2)  .  All  of  the  ( k  —  2)  remaining  main 
effects  must  be  sent  to  pins  or  minus  the  other  (k  —  2)  viable  2-factor  interactions. 
That  is,  there  are  (2k)(2(k  —  1))(2 k~2(k  —  2)!)  =  2 kk\  distinct  possibilities. 

(Case  3:  2-factor  interaction,  main  effect)  Suppose  g(xi)  =  ix^j  for  i  <  j  and 
i,j  E  {l,...,k}  and  g(x2)  =  ±x^  for  l  e  {i,j}  ■  Then  there  are  2(^)  possibilities 
for  gfei)  and  2(2)  for  (?(x2)  .  All  of  the  ( k  —  2)  remaining  main  effects  must  be  sent 
to  plus  or  minus  the  other  (k  —  2)  viable  2-factor  interactions.  That  is,  there  are 
(2g))(2(2))(2 k~2(k  -  2)!)  =  2k(k)(2 )(k  -  2)!  distinct  possibilities. 

(Case  4:  2-factor  interaction,  2-factor  interaction)  Suppose  g(xi)  =  ±x*j  for  i  <  j 
and  i,j  E  {1, . . .  ,/c}  and  g(x2)  =  ±x,  for  l  E  {(l,i),  (1,  j),  •  •  •  >  (®  — 1, i),  (i-l,j),  (i,i  + 
1),  (i+l,j),..  • ,  (i,j- 1),  (j-1,  j),  (i,j+ 1),  (j,j+ !),••  • ,  (i,k),  ( j,k )}  .  Then  there  are 
2Q)  possibilities  for  ^(xi)  and  2(2k  —  4)  for  ^(x2)  .  All  of  the  (k  —  2)  remaining  main 
effects  must  be  sent  to  plus  or  minus  the  other  (k  —  3)  viable  2-factor  interactions  and 
the  lone  viable  main  effect  xpjjn;  .  That  is,  there  are  (2(k))(2(2k  —  A))(2k~2(k  —  2)!)  = 
2k(k)(2 k  —  4)(k  —  2)!  distinct  possibilities. 
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Therefore,  the  total  number  of  possibilities  for  all  cases  is 


2‘(*W-4)(*.-2)! 

=  2t(k-2)\(2  *(*-l)+Q +  (*)(*: -2)  ) 

=  2fc(fe  -  2)!  (^2  -  1)  +  Q)  (*  -  !)] ) 

=  2*(fc- l)l(2  *+Q  ) 

=  2k(k  —  l)\(2k  +  k(k  —  1)) 

=  2k{k)\(2  +  k-  1) 

=  2k{k  +  1)! 

Thus,  g  is  one  of  at  most  2 k{k  +  1)!  elements  in  G  . 

Theorem  3  Let  k  >  4  .  Then 


2kk\  +  2kk\  +  2k  (  )  {2){k  —  2)!  + 


G  =  Sk  x  Sk+1 


Proof  Let  R  =  (pi, . . . ,  pk)  where  pi  acts  on  the  full  factorial  design  by  sending 

(xi, . . . ,  Xj, . . . ,  Xfc)  to  (xx  j, . . . ,  Xj, . . . ,  Xj  fc)  for  i  =  1, _ ,  k  .  Note  that  elements  of 

R  preserve  the  full  factorial  design  as  well  as  R,ow(M)  ,  so  R  <  G  .  Furthermore,  for 
i  =  1, . . . ,  k  ,  p'1  =  Pi  ■  For  any  distinct  i ,  j  G  {1, . . . ,  k}  ,  PiPjPi  simply  permutes 
Xj  and  Xj  within  the  full  factorial  design,  so  clearly  Sk  <  R  .  Now  we  see  that 
PjPiSk  =  PiSk  for  distinct  i,  j  G  {1, . . . ,  k}  ,  so  there  are  exactly  k  +  1  left  cosets  of  Sk 
within  R  ,  and  together  these  constitute  the  entirety  of  R  .  Hence,  R  =  Sk+1  .  Letting 
<f)  G  Sk  be  arbitrary,  we  note  that  for  any  i  —  1, . . . ,  k  ,  pj1  =  p*  ,  p~x^>pi  =  <ft'  where 
(//  G  Sk  ,  and  p~X(\)' pi  =  (p  .  Together  with  this  information,  Lemma  4  makes  it  clear 
that  Sk  <  S%  x  Sk+ 1  .  Now  S'g  x  <  G  ,  and  \Sk  x  =  \Sk \  l^fe+il  =  2k(k+  1)! 
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which  is  the  upper  bound  for  (S'  ,  so  G  must  be  exactly  5^  x  Sk+ 1  • 
Corollary  2  Let  k  >  4  .  Then 


\G\  =  2k(k  +  1)! 

Proof  As  a  direct  result  of  Theorem  3,  we  have  (S'  =  |5f  x  5fc+i  =  2k{k  +  1)!  . 

Note  that  when  k  =  3  ,  the  result  above  does  not  hold.  For  example,  consider  the 
permutation  g  e  G  such  that 

g(-xi)  =  0.5x1i2  +  0.5x1j3  +  0.5x2  -  0.5x3 
5»(x2)  =  0.5x1j2  +  0.5x1i3  -  0.5x2  +  0.5x3 
.r/(x3)  =  Xi 

Because  this  permutation  sends  main  effects  to  forms  other  than  those  which  were 
viable  for  k  >  4  ,  we  conclude  (S'  >  23(3+l)!  =  192  .  This  observation  is  corroborated 
by  the  Geyer  [2]  algorithm  and  GAP  [1],  which  prove  that  in  this  case  |(S|  =  1152  , 
and  G  =  (54  x  54)  x  52  . 
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IV.  TOPICS  FOR  FUTURE  RESEARCH 


It  is  easy  to  see  that  S2 1  Sk  is  always  a  subgroup  of  G  ,  and  Sf  x  Sk+i  is  a  subgroup 
when  t  is  even.  Perhaps  the  entirety  of  G  can  be  found  for  OA (N,  k,  2,  t)  with  t  >  2  . 
Also,  finding  the  LP  relaxation  symmetry  group  of  the  ILP  formulation  in  [3]  of  an 
OA (N,  k,  s,  t )  for  s  >  2  is  an  open  problem.  In  this  case,  it  is  easy  to  see  that  Ss  l  Sk 
is  always  a  subgroup  of  this  group. 
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Appendix  A.  R  CODE  AND  CASES 


################################################################################ 

#  Function-  expandcases  -  generates  linear  combinations  with,  positive  1st  term 

#  Input (s)-  1st  -  binary  representation  of  main  effects  and  interactions 

#  Output (s)-  newlst  -  linear  combinations  with  positive  1st  term 
################################################################################ 

expandcases<-function(lst) { 
ll<-length(lst) 
newlst<-list () 
r<-l 

aa<-as . matrix (expand. grid ( c ( 1 , -1) ,c(l,-l) ,c(l,-l))) 
aaOcbind  ( 1 ,  aa) 
for  (j  in  1:11){ 
for  (i  in  1:8){ 
aa2<-rbind(aa[i,] ,1st  [[j]] ) 
dimnames (aa2) [ [2] ] <-NULL 
newlst [r] <-list (aa2) 
r<-r+l 
» 

return (newlst) 

> 

################################################################################ 

#  Function-  allcheckbinvector  -  checks  viability  of  linear  combinations 

#  Input (s)-  1st  -  all  linear  combinations  to  be  checked 

#  Output (s)-  displays  Is  for  viable  combinations  and  Os  otherwise 
################################################################################ 

allcheckbinvector<-function(lst) { 

lc-length(lst) 

for  (i  in  1:1){ 

print (checkbinaryvector (1st [ [i]  ]  ) ) 

> 

> 

################################################################################ 

#  Function-  checkbinaryvector  -  checks  viability  of  a  linear  combination 

#  Input (s)-  newcik  -  the  linear  combination  to  be  checked 

#  Output (s)-  1  if  combination  is  viable,  0  otherwise 
################################################################################ 

checkbinary vector<-funct ion (newcik) { 

newcikmat<-newcik [— 1 ,  ] 

newcikcof f <-newcik [1 ,  ] 

pp<-dim(newcikmat) 

cols<-pp [2] 

pp<-pp[l] 

full<-c ("expand.gr id (c(l , -1) ") 
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for  (j  in  1 : (pp-1) ){ 
full<-paste (full , " ,c(l,-l) ") 

> 

f ull<-paste (full , " ) " ) 
full<-parse (text  =  full) 
full<-eval (full) 

temp<-genprod(full ,newcikmat [ , 1]  ) 
for  (i  in  2:cols){ 

temp<-cbind (temp, genprod (full ,newcikmat [, i] ) ) 

} 

f invec<-as .matrix(temp [ , 1] *newcikcoff [1] ) 
for  (i  in  2:cols){ 

f invecC-f invec+newcikcof f [i] *temp [ , i] 

> 

f invecC-as .matrix(as . integer (finvec/2) ) 

if (min(f invec)==-l  &  max(finvec)==l){ 
return(l)}  else-Creturn(O)} 

} 

################################################################################ 

#  Function-  genprod  -  computes  Hadamard  product 

#  Input (s)-  full  -  2-level  (+/-1)  full  factorial  design 

#  ttl  -  indicator  of  which  main  effects  are  to  be  multiplied 

#  Output(s)-  as .matrix(outfullcheck)  -  +/-1  form  of  main  effect  or  interaction 
################################################################################ 

genprodofunct  ion  (full ,  tt  1)  { 
outfullcheckofull  [,  1]  ~  (ttl  [1]  ) 
pp<-length (ttl) 
for  (j  in  1 : (pp-1) ){ 

outfullcheck<-outfullcheck*full  [, (j+1)] ~ (ttl  [(j+1)] ) 

} 

return (as . matrix (out full check) ) 

> 

################################################################################ 

#  Cases  with  1  two-factor  interaction 

################################################################################ 

al2.1.2.3<-cbind(c(l,l,0) ,c (1,0,0) ,c (0,1,0) ,c (0,0,1)) 
cl2 . 1 . 2 . 3<-expandcases (list (al2 .1.2.3)) 
allcheckbinvector (cl2 .1.2.3) 

al2 . 1 .3 .4<-cbind(c(l , 1,0,0) ,c(l, 0,0,0) ,c (0,0, 1,0) ,c (0,0, 0,1)) 
cl2 . 1 . 3 . 4<-expandcases (list (al2 .1.3.4)) 
allcheckbinvector (cl2 .1.3.4) 


24 


al2.3.4.5<-cbind(c(l,l,0,0,0) ,c (0,0, 1,0,0) ,c (0,0, 0,1,0) ,c (0,0, 0,0,1)) 
cl2 . 3 . 4 . 5<-expandcases (list (al2 .3.4.5)) 
allcheckbinvector (cl2 .3.4.5) 

################################################################################ 
#  Cases  with  2  two-factor  interactions 

################################################################################ 

al2 . 13<-cbind(c(l ,1,0, 0,0,0) ,c(l, 0,1, 0,0,0)) 
al2 . 34<-cbind(c (1 , l,0,0,0,0),c(0J0,l,l,0,0)) 

al2 . 13 . l<-cbind(al2 . 13, c(l, 0,0, 0,0,0)) 
al2 . 13 . 2<-cbind(al2 . 13, c(0, 1,0, 0,0,0)) 
al2 . 13 . 4<-cbind(al2 . 13, c(0, 0,0, 1,0,0)) 
al2.34. l<-cbind(al2.34,c(l,0,0,0,0,0)) 
al2 . 34 . 5<-cbind(al2 . 34 , c (0 , 0 , 0 , 0 , 1 , 0) ) 

al2 .13.1. 2<-cbind(al2 . 13.1, c(0, 1,0, 0,0,0)) 
al2 .13.1. 2<-al2 . 13. 1.2  [c (1:3),] 
cl2 . 13 . 1 . 2<-expandcases (list (al2 .13.1.2)) 
allcheckbinvector (cl2 .13.1.2) 

al2 . 13 . 1 ,4<-cbind(al2 . 13. l,c(0, 0,0, 1,0,0)) 
al2 .13.1 ,4<-al2 . 13 . 1 .4 [c(l :4) ,] 
cl2 . 13 . 1 . 4<-expandcases (list (al2.13.1.4)) 
allcheckbinvector (cl2 .13.1.4) 

al2 .13.2. 3<-cbind(al2 . 13.2, c(0, 0,1, 0,0,0)) 
al2 . 13 . 2 . 3<-al2 . 13. 2. 3 [c (1:3) ,] 
cl2 . 13 . 2 . 3<-expandcases (list (al2 .13.2.3)) 
allcheckbinvector (cl2 . 13.2.3) 

al2 . 13 . 2 .4<-cbind(al2 . 13.2, c(0, 0,0, 1,0,0)) 
al2 . 13 . 2 ,4<-al2 . 13 . 2 . 4 [c (1 : 4) ,  ] 
cl2 . 13 . 2 ,4<-expandcases (list(al2.13.2.4)) 
allcheckbinvector (cl2 . 13.2.4) 

al2 . 13 .4. 5<-cbind(al2 . 13.4, c(0, 0,0, 0,1,0)) 
al2 . 13 .4. 5<-al2 . 13. 4. 5 [c (1:5),] 
cl2 . 13 . 4 . 5<-expandcases (list (al2 .13.4.5)) 
allcheckbinvector (cl2 .13.4.5) 

al2 .34.1. 2<-cbind(al2 . 34. 1 , c(0, 1 ,0, 0,0,0) ) 
al2 .34.1. 2<-al2 . 34 . 1 . 2 [c (1 : 4) , ] 
cl2 . 34 . 1 . 2<-expandcases (list (al2 .34.1.2)) 
allcheckbinvector (cl2 .34.1.2) 

al2 .34.1. 3<-cbind(al2 . 34 . 1 , c(0,0, 1 , 0,0,0)  ) 
al2 .34.1. 3<-al2 . 34 . 1 . 3 [c (1 : 4) ,  ] 
cl2 . 34 . 1 . 3<-expandcases (list (al2 .34.1.3)) 
allcheckbinvector (cl2 .34.1.3) 
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al2 . 34 . 1 . 5<-cbind(al2 . 34.1, c(0, 0,0, 0,1,0)) 
al2 .34.1. 5<-al2 . 34.1. 5 [c (1:5) ,] 
cl2 . 34 . 1 . 5<-expandcases (list (al2 .34.1.5)) 
allcheckbinvector (cl2 .34.1.5) 

al2 . 34 . 5 . 6<-cbind(al2 .34.5,0(0,0,0,0,0,1)) 
al2 . 34 . 5 . 6<-al2 . 34. 5. 6  [c (1:6)  ,] 
cl2 . 34 . 5 . 6<-expandcases (list (al2 .34.5.6)) 
allcheckbinvector (cl2 .34.5.6) 

################################################################################ 
#  Cases  with  3  two-factor  interactions 

################################################################################ 

al2.13<-cbind(c (1,1, 0,0, 0,0,0) ,c (1,0, 1,0, 0,0,0)) 
al2.34<-cbind(c (1,1, 0,0, 0,0,0) ,c (0,0, 1,1, 0,0,0)) 

al2 . 13 . 14<-cbind(al2 . 13,c(l,0,0,l,0,0,0)) 
al2 . 13 . 23<-cbind(al2 . 13, c(0 , 1 , 1 ,0,0,0 ,0) ) 
al2 . 13 . 24<-cbind(al2 . 13, c(0 , 1 ,0, 1 ,0,0 ,0) ) 
al2.13.45<-cbind(al2. 13, c (0,0, 0,1, 1,0,0)) 
al2 . 34 . 56<-cbind(al2 . 34, c (0 , 0 , 0 , 0 , 1 , 1 , 0) ) 

al2 . 13 . 14. l<-cbind(al2 . 13.14, c(l, 0,0, 0,0, 0,0)) 
al2 . 13 . 14. l<-al2 . 13 . 14. 1 [c(l : 4) ,] 
cl2 . 13 . 14. l<-expandcases (list (al2 .13.14.1)) 
allcheckbinvector (cl2 .13.14.1) 

al2 . 13 . 14. 2<-cbind(al2 . 13. 14, c (0,1, 0,0, 0,0,0)) 
al2 . 13 . 14. 2<-al2 . 13. 14. 2 [c (1:4) ,] 
cl2 . 13 . 14. 2<-expandcases (list (al2 .13.14.2)) 
allcheckbinvector (cl2 . 13. 14.2) 

al2 . 13 . 14. 5<-cbind(al2 .13. 14, c(0, 0,0, 0,1, 0,0)) 
al2 . 13 . 14. 5<-al2 . 13 . 14. 5 [c(l : 5) ,] 
cl2 . 13 . 14. 5<-expandcases (list (al2 .13.14.5)) 
allcheckbinvector (cl2 .13.14.5) 

al2 .13.23. l<-cbind(al2 . 13.23, c(l, 0,0, 0,0, 0,0)) 
al2 .13.23. l<-al2 . 13 . 23. 1 [c(l : 3) ,] 
cl2 . 13 . 23 . l<-expandcases (list (al2 .13.23.1)) 
allcheckbinvector (cl2 .13.23.1) 

al2 .13.23. 4<-cbind(al2 . 13 . 23 , c (0 , 0 , 0 , 1 , 0 , 0 , 0) ) 
al2 .13.23. 4<-al2 . 13 . 23 . 4  [c ( 1 : 4) , ] 
cl2 . 13 . 23 . 4<-expandcases (list (al2 .13.23.4)) 
allcheckbinvector (cl2 .13.23.4) 

al2 . 13 . 24. l<-cbind(al2 .13. 24, c(l, 0,0, 0,0, 0,0)) 
al2 .13.24. l<-al2 . 13 . 24. 1 [c(l :4) ,] 
cl2 . 13 . 24. l<-expandcases (list (al2 .13.24.1)) 
allcheckbinvector (cl2 .13.24.1) 
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al2 . 13 . 24. 3<-cbind(al2 . 13. 24, c(0, 0,1, 0,0, 0,0)) 
al2 .13.24. 3<-al2 . 13 . 24 . 3 [c ( 1 : 4) , ] 
cl2 . 13 . 24. 3<-expandcases (list (al2 .13.24.3)) 
allcheckbinvector (cl2 . 13.24.3) 

al2 . 13 . 24. 5<-cbind(al2 . 13. 24, c (0,0, 0,0, 1,0,0)) 
al2 .13.24. 5<-al2 . 13 . 24 . 5 [c ( 1 : 5) , ] 
cl2 . 13 . 24. 5<-expandcases (list (al2 .13.24.5)) 
allcheckbinvector (cl2 . 13.24.5) 

al2 . 13 . 45 . l<-cbind(al2 . 13. 45, c(l, 0,0, 0,0, 0,0)) 
al2 . 13 .45 . l<-al2 . 13 .45 . 1  [c(l : 5) ,] 
cl2 . 13 . 45 . l<-expandcases (list (al2 .13.45.1)) 
allcheckbinvector (cl2 .13.45.1) 

al2 .13.45. 2<-cbind(al2 . 13 .45, c(0, 1 , 0,0,0,0 ,0) ) 
al2 . 13 . 45 . 2<-al2 . 13 . 45 . 2 [c (1 : 5) ,] 
cl2 . 13 . 45 . 2<-expandcases (list (al2 .13.45.2)) 
allcheckbinvector (cl2 . 13.45.2) 

al2 .13.45. 5<-cbind(al2 . 13 .45, c(0 ,0, 0,0, 1 ,0 ,0) ) 
al2 .13.45. 5<-al2 . 13 . 45 . 5  [c ( 1 : 5) , ] 
cl2 . 13 . 45 . 5<-expandcases (list (al2 .13.45.5)) 
allcheckbinvector (cl2 . 13.45.5) 

al2 .13.45. 6<-cbind(al2 . 13 . 45, c(0 ,0, 0,0,0, 1 ,0) ) 
al2 . 13 . 45 . 6<-al2 . 13. 45. 6 [c (1:6) ,] 
cl2 . 13 . 45 . 6<-expandcases (list (al2 .13.45.6)) 
allcheckbinvector (cl2 . 13.45.6) 

al2 .34.56. l<-cbind(al2 .34.56,0(1,0,0,0,0,0,0)) 
al2 .34.56. l<-al2 . 34. 56 . 1 [c(l : 6) ,] 
cl2 . 34 . 56 . l<-expandcases (list (al2 .34.56.1)) 
allcheckbinvector (cl2 .34.56.1) 

al2 .34.56. 7<-cbind(al2 . 34 . 56 , c (0 , 0 , 0 , 0 , 0 , 0 , 1) ) 
cl2 .34.56. 7<-expandcases (list (al2 .34.56.7)) 
allcheckbinvector (cl2 .34.56.7) 

################################################################################ 
#  Cases  with  4  two-factor  interactions 

################################################################################ 

al2 . 13<-cbind(c (1 ,l,0,0,0,0,0,0),c(l, 0,1, 0,0, 0,0,0)) 
al2 . 34<-cbind(c (1 ,l,0,0,0,0,0,0),c(0, 0,1, 1,0, 0,0,0)) 

al2 . 13 . 14<-cbind(al2 . 13, c(l, 0,0, 1,0, 0,0,0)) 
al2 . 13 . 23<-cbind(al2 . 13, c(0, 1,1, 0,0, 0,0,0)) 
al2 . 13 . 24<-cbind(al2 . 13, c(0, 1,0, 1,0, 0,0,0)) 
al2 . 13 . 45<-cbind(al2 . 13, c(0, 0,0, 1,1, 0,0,0)) 
al2 . 34 . 56<-cbind(al2 . 34, c(0, 0,0, 0,1, 1,0,0)) 
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al2 .13.14. 15<-cbind(al2 . 13. 14, c(l, 0,0, 0,1, 0,0,0)) 
al2. 13. 14. 15<-al2. 13. 14. 15  [c(l :5) ,] 
cl2 . 13 . 14. 15<-expandcases (list (al2 . 13 . 14 . 15) ) 
allcheckbinvector (cl2 .13.14.15) 

al2 . 13 . 14. 23<-cbind(al2 . 13. 14, c (0,1, 1,0, 0,0, 0,0)) 
al2 .13.14. 23<-al2 . 13 . 14 . 23 [c (1 : 4) , ] 
cl2 . 13 . 14. 23<-expandcases (list (al2 . 13 . 14 . 23) ) 
allcheckbinvector (cl2 . 13. 14.23) 

al2 . 13 . 14. 25<-cbind(al2 . 13. 14, c (0,1, 0,0, 1,0, 0,0)) 
al2 .13.14. 25<-al2 . 13 . 14 . 25 [c (1 : 5) ,  ] 
cl2 . 13 . 14. 25<-expandcases (list (al2 . 13 . 14 . 25) ) 
allcheckbinvector (cl2 . 13. 14.25) 

al2.13.14.56<-cbind(al2.13.14,c(0,0,0,0,l,l,0,0)) 
al2 .13.14. 56<-al2 . 13 . 14 . 56  [c (1 : 6) , ] 
cl2 . 13 . 14. 56<-expandcases (list (al2 . 13 . 14 . 56) ) 
allcheckbinvector (cl2 . 13. 14.56) 

al2.13.23.45<-cbind(al2. 13.23,c(0,0,0,l,l,0,0,0)) 
al2 . 13 . 23 . 45<-al2 . 13 . 23 . 45  [c(l : 5) , ] 
cl2 . 13 . 23 . 45<-expandcases (list (al2 .13.23. 45) ) 
allcheckbinvector (cl2 .13.23 .45) 

al2 . 13 . 24. 35<-cbind(al2 . 13. 24, c (0,0, 1,0, 1,0, 0,0)) 
al2 .13.24. 35<-al2 . 13 . 24 . 35 [c (1 : 5) ,  ] 
cl2 . 13 . 24. 35<-expandcases (list (al2 . 13 . 24 . 35) ) 
allcheckbinvector (cl2 . 13.24.35) 

al2 . 13 . 24. 56<-cbind(al2 . 13. 24, c (0,0, 0,0, 1,1, 0,0)) 
al2 . 13 . 24. 56<-al2 . 13 . 24. 56  [c(l : 6) ,] 
cl2 . 13 . 24. 56<-expandcases (list (al2 . 13 . 24 . 56) ) 
allcheckbinvector (cl2 . 13.24.56) 

al2.13.45.56<-cbind(al2. 13.45,c(0,0,0,0,l,l,0,0)) 
al2 .13.45. 56<-al2 . 13 . 45 . 56  [c ( 1 : 6) , ] 
cl2 . 13 . 45 . 56<-expandcases (list (al2 . 13 .45 . 56) ) 
allcheckbinvector (cl2 . 13.45.56) 

al2 .13.45. 67<-cbind(al2 .13.45,0(0,0,0,0,0,1,1,0)) 
al2 . 13 .45 . 67<-al2 . 13 .45 . 67 [c(l :7) ,] 
cl2 . 13 . 45 . 67<-expandcases (list (al2 . 13 .45 . 67) ) 
allcheckbinvector (cl2 . 13.45.67) 

al2 . 34 . 56 . 78<-cbind (al2 . 34. 56, c(0, 0,0, 0,0, 0,1,1)) 
cl2 . 34 . 56 . 78<-expandcases (list (al2 .34.56. 78) ) 
allcheckbinvector (cl2 .34.56. 78) 
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